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Nonequlibrium phase transition of an open Takayama-Lin Liu-Maki chain coupled with 
two reservoirs is investigated. We will show that solitons connecting two uniform phases are 
possible, and the amplitude of solitons obeys the same self-consistent equation as that of 
uniform phases. 

§1. Introduction 

The study of nonequihbrium transport of electrons is essential to understanding 
the underlying physics of nano devices. Because of the development of nano devices, 
theoretical description of quantum transport may lead to experimentation to test 
the validity of quantum theories. As such, the pure science of quantum transport 
holds significant value in contemporary mesoscopic physics. 

A common approach to the study of quantum transport in low-dimensional sys- 
tems is putting systems with infinitely extended reservoirs into contact with one 
another. However, an analytical approach to such systems is usually difficult, and 
the majority of works are performed numerically. Nevertheless, recent progress in 
analytical frameworks allows us to study several systems at nonequilibriumn steady 
state (NESS). The algebraic approach is promising for the study of NESS in such 
systems,'3S''3 and NESS is studied in a small number of cases, e.g., suppression of 
Fano-Kondo plateau in Aharonov-Bohm rings,® nonlinear conductance of a solvable 
model of the Kondo effect j'S' long range correlations in the XY model,!^''^ and a new 
phase transition of the Takayama-Lin Liu-Maki (TLM) model.l^ 

Using a result of the algebraic approachP''SJ''I^''n'' (Eq. ()3Tp in §3), we studied 
NESS of the TLM modeP as a representative example of a 1-D system exhibit- 
ing collective orders. In that work, we discussed isotropic phases (i.e., uniformly 
dimerized chains) and the suppression of order under the application of current and 
bias voltage. Also, we found orders to be multi-valued functions of bias voltage and 
single- valued functions of current. Since the continuity equation implies that current 
is spatially uniform, the existence of a solution which connects two uniform domains 
with the equal amplitude is suggested. 

Indeed, it is known that only solitons (amplitude kinks) connecting two domains 
with equal amplitudes are possible at equilibrium 

In this paper, we discuss such solitons at NESS. As in the equilibrium case, the 
corresponding fermionic spectrum consists of continuum states and a midgap state 
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at the center of the energy gap. 

However, since our system is connected to reservoirs, reservoir fermions inside 
the energy gap should be considered. We will show that in spite of the existence 
of the midgap state, reservoir fermions carrying the wave number inside the energy 
gap do not contribute to the gap equation, resulting in the same gap equation as the 
isotropic phases. 

§2. Open TLM Model and the Equation of Motion 

The TLM modeP^ is a continuum limit of a tight-binding model (the SSH lat- 
tice) for polyacetylene proposed by Su, Schrieffer, and Heegei'^''^ which describes 
charge density wave commensurate with the lattice. 

The TLM chain consists of two fermionic fields, d{x) and e{x), and the quantized 
local lattice distortion A{x). Our Hamiltonian is given by (see appendix A of ref. [8]) 
for its derivation) 



Hs= dx^\x) 
Jo 



-ihvay— — h A(x)ax 
ox 



•P'(x) 



1 f . 
H r / dx 

ZTTnvX Jq 



A{xf + \n{xf 



where ^{x) = {d{x) , e{x))'^ is the two-component fermionic field, i is the length 
of the system, v is the Fermi velocity, ax and are the Pauli matrices, A is the 
dimensionless coupling constant, ujq is the phonon frequency and 77(x) corresponds to 
the momentum conjugate to A{x). Nonvanishing equal-time commutation relations 
among those operators are 

{d{x),d{yy} = {e{x), e(y)t} = 5{x - y) , \A{x), i7(y)] = ih^TT\vul5{x - y) , 

where {A, B} = AB + BA and [^4, B] = AB — BA. As the system is finite, fermionic 
waves are reflected back at the edges. This effect is taken into account by a boundary 
condition: 

d(0) = , e{C) = . 
The reservoirs are described by 

= j dk{hujkLa\i^akL + ^kRO^kR^'kR} i 

where afci,, (z/ = L,i?) stands for the annihilation operators of reservoir fermions 
with momentum fc, hw^u represents their energies, and their nonvanishing anticom- 
mutation relations are written {afci/,a|,/^} = 5{k — k'). The TLM chain-reservoir 
interaction is assumed to be 

V = j dk |^ffce"''(0)afcL + hwkd\t)akR + [h.c] 
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where Vk and Wk stand for the couphng matrix elements. Then, the Hamiltonian of 
the whole system is given by 

H = Hs + V + Hb . (2-1) 
From (|2-1|) . the lattice distortion A is found to obey the following equation of motion: 
d^A{x,t) _ dn{x,t) 



dt 



-UJn 



'o ( A{x, t) + TThvX'F^x, t)ax'^{x, t) 



2-2 



38. NESS Mean-Field Approximation 



In this section, we discuss mean-field NESS of the TLM model. NESS under the 
mean-field approximation (• • • )oo can be characterized as a state satisfying Wick's 
theorem with respect to the incoming fields q^^ of a^u for the mean-field Hamiltonian 



Hf^ + V + HB 



H 



MF 
5 



d 



-ihvay— — h A(x)ax 
ox 



'Fix) 



Namely, they are defined as the solution of 
1 



h 



-00 



(3-1) 



If asymptotic fields are complete, the incoming fields satisfy: 

(«L"fcv)oo = fu{fu^ku)S{k - k') , 

where fiuj^u is the single-particle energy of wave number k, fu{x) = \j{e^^~^^^l'^^ + 1) 
is the Fermi distribution function, Ty is the initial temperature, and [ly is the initial 
chemical potential of the reservoirs [y = L,R). 

The self-consistent equation for the lattice distortion A{x) = {A{x))oo is derived 
from the equation of motion (|2-2p for the lattice distortion and the time-independence 
of A{x). Then, 1^ leads to 

-ld'^A{x) 







UJn 



A{x) +7rhvX{'l'\x)ax^{x))c 



(3-2) 



Hereafter, we study a case where the lattice distortion is a soliton (amplitude kink) 
connecting two domains with the equal amplitude: A{x) = Aotanh Ks{x — a). More- 
over, we restrict ourselves to the case of Aq = hvKg, which corresponds to cases 
in which creation energy of a soliton is minimum, and in which only such solitons 
satisfy the self-consistent equation at equilibrium(see ref. [T2|) ). In this case, the 
self-consistent equation (j3-2p reads as 



/ dcul 
J —00 I 



T ? f .h{x;uj)^axh{x;uj) 

\A^{uj)\^ ^^(^^ 

h{x;uj)^axh{x;uj) 
+Im?7_(w) 1^ ^^^p 



A{x) 
ttHvX 



(3-3) 
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where the integral is taken over by the reservoir fermions' energy, A-{uj), h{x;uj) and 
h{x; oj) are auxihary functions defined by (the green function is given in appendix Rl 
and relations between the auxiliary functions and the green function was shown in 
appendix B of ref. |8|)) 



DA_{uj) = (hvK) + {A{1) - Z\(0)) - hvKi.T]. 



UJ 



COS k£ 



+ 



hi{x) 
h2{x) 



hi{x) 
h2{x) 



1 



1 



z\(o) + + ^-) + 

cos n{x — I) 
1 ( A{x) I 



7?„ + A{(.)^^ri^ 



sin kI 



-2^{A{x)-A{0)} + hvKi 
—hvK 



COS KX 



- dk' — 

V J UJ 



V±{^) = - 

iOk'L ± V 




hio + A{e)r]+ 



sin k{x — i) 



{A{0) + hio)i 
hnv cos kI + A{0) sin kI 



sm KX 



Equation (j3-3p is the self-consistent equation for the order parameter A{x). Note 
that we use a convention Im^_(w) = {lmjj-{Lj) = 0) for u outside the range of tOku- 
We further note that in accordance with Aq = hvKs, a^L depends on position x and 
a^R does not depend on position in case of Aq = hvKs. 

§4. Gap Equation for Solitons 

4.1. Self- consistent Equation 

In this section, we show that the left-hand side of Eq. (l3-3p is proportional to 
A(x) = Z\otanhK;s(x — a) for a long chain, and derive the self-consistent equation 
for solitons' amplitude Aq. 

Hereafter, we consider a case where the zero-bias chemical potential is located 
at the band center of the TLM chain. In other words, the chain is half-filled. This 
case corresponds to fkoku = (/i|fc|)^/(2m^) — eo^ {ly = L,R) with and being, 
respectively, the effective mass and the zero-bias chemical potential measured from 
the band bottom of the reservoirs. In order to prevent the increase of electrostatic 
energy, the chemical potentials of the reservoirs should be chosen so that = 
— /i/j = —eV/2 with V signifying bias voltage, and e the elementary charge (for detail, 
see appendix C of ref. Ej)). As is well known,!!^ the energy cutoff fkUc is necessary for 
the TLM model and we assume that cqi, — e\V\/2 > hwc- The integration interval 
of left-hand side of Eq. ()3-3p should be replaced by (— Wc^c)- In contrast to the 
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isotropic dimerization, integrant of |?ia;| < Aq should be carefully treated due to the 
existence of the midgap state. In the following two subsections, we evaluate left-hand 
side, of Eq. (l3^ . 



4.2. Reservoir Fermions with Energy inside Energy Gap 

In this subsection, the left-hand side of Eq. ()3-3p is evaluated in the energy 
range of I/k^I < Aq, which corresponds to the energy gap of the system. This 
appears exponentially small, as in the case of uniform phases. However, there exists 
a midgap state at the center of the gap, and the exponentially large term with respect 
to i in the denominator DA^ becomes comparable to the exponentially small term 
with respect to £. Therefore, we employ the new integrating variable ey = hw with 
e = exp(— Ks^) to scale up the behavior near w = 0. Then, DA- is evaluated by 



2D A- = - 

y 



2 



+ 0{e) 



where the first three terms are exponentially large, and the last term is exponentially 
small in the original integrating variable lo. Since we are considering a long chain, 
we focus on the center of the chain. For this purpose, we further introduce a new 
variable x = 1/2 + 5x. It is thus easy to prove that after rewriting the left-hand side 
of Eq. (|3-3p with y and 5x, terms on the order of and e'^ vanish. As a result, we 
conclude that the contribution of reservoir fermions with energy inside the energy 
gap to the gap equation is exponentially small, and thus is negligible. Note that 
this is in contrast to the result of polarons which only appear out-of-equilibrium due 
to the contribution of reservoir fermions carrying energy inside the gap to the gap 
equation. In that case, reservoir fermions near the midgap contribute to the gap 
equation and induce new excitation to the spinless TLM model (this aspect will be 
discussed elsewhere). 

4.3. Reservoir Fermions Carrying Energy inside Continuum State 

In this subsection, the left-hand side of Eq. (j3-3p is evaluated in the energy range 
of > Aq. Since the dominant terms of the left-hand side of Eq. (j3-3p do not 
change regardless of energy, and because we consider a case where the center of the 
soliton is far from the chain end, i.e., a = 0{i), {I — a) = 0{i), the lattice distortion 
at the chain ends, i.e., A{0) and A{1), are approximated by —A and A. Then, e.g., 
the second term of the left-hand side of Eq. ()3-3p reads as 



f h{x;uj)^a^h{x;uj) 
/ d^ wlmr?_(w) /ij(^) 

JAo/h<H<u;, 

duJ — 1 ^ . , , — — r77r\ ^2\^] sm Ikx + C/2(w) cos Ikx 
Ao/h<\u^\<uj, v\D{u:)A-{uj)Y I 

+ hA{x)(u: - ^Re^lH (4-1) 

where C2{u!) and S2iio) are defined by 

C2(cj) = 2 {AoA{x) - [hvnf] ReC-(w) - nujA{x)\i.{uj)\^ 
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2Ao{hvKf + A{x) {{hvKf - Al} 



+ 



2 



5'2(u;) = 2{hvK) {A{x) + Aq) Re^_(a;) - huj{hvK)\£,-{uj 

TVjJ 

Moreover, by applying the Riemann-Lebesgue lemma, (|4-ip is found to be 

f , h[x]u])^axh{x]U}) 

JAo/h<\i,j\<UJa l^-l'^jl 



where corresponds to a dominant term of the Fourier coefficient of l/|L'yl_(tij)|^ 

and is given by 



Co(w) = h'^VK 



2A 

Im,^_(w) ^uj -\ —Kerj^{LL!) + 6<j|r/_(L<j)| 



2A ~ ~ 
+Imr7„(Lt;) \ ijJ ^Re,^_(a;) + lij|^_(l<j)| 



By a similar argument, when the TLM chain is long enough, the first term of the 
left-hand side of Eq. (j3-3p is proportional to A{x), and Eq. (j3-3p leads to 



A 

where 



Z\o = or (4-2) 
S{A^,V,Tl,Tr) , (4-3) 



^2 = y d^|^-^(a; + ^Re^_H+^|r/-MI j/L(M 



UinI I I 



^ Imr?_(a;) / 2Aq 



Co(u;) V ?i 



Ree_H+c^|e-Mp)/R(M 



When the TLM chain couples with two reservoirs at temperature T = Tl = Tr, an 
assumption = — ^_(u;) simplifies the self-consistent equation (j4-3p : 

1 f'"- sinh(e/r) 

A J|zio| - Z\o cosh(l^) + cosh(e/r) 

We would like to remark that the amplitude of solitons Aq and isotropic phases obey 
the same self-consistent equation (Eq. (j4-4p ). though another assumption rj-{u;) = 
^_(a;) is required for isotropic phases (this difference is induced by the phase shifts 
of the solitons). 
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4.4. Fermionic Current 



The fermionic current at x in the TLM chain is given by 

J{x) = -ev>F^{x)ay>F{x) . 



Since current is carried by fermions with energies near the fermi level, we approximate 
(,-{uj), byi^_(u;) = rj-{io) = iF/v. By applying the Riemann-Lebesgue lemma, 

its NESS average is reduced to 



where Go = e^vF / {-nhiv"^ is the conductance in the normal phase. We would 

like to note that the effect of fermions carrying energy inside the energy gap is 
exponentially small with respect to the chain size ^, and is negligible. 



We have studied nonequilibrium Peierls transition in the Takayama-Lin Liu-Maki 
chain connected to two reservoirs at different chemical potentials (their difference cor- 
responds to bias voltage) by combining a mean-field approximation and the formula 
(j3-ip . which is an outcome of the algebraic field approach to nonequilibrium statis- 
tical mechanics. We show that amplitude of solitons and isotropic lattice distortion 
obey the same self-consistent equation ()4-4p after certain approximations; however, 
solitons require another condition of tunneling between systems and reservoirs be- 
cause of the effect of phase shifts. We only consider solitons whose creation energy 
is minimum {Aq = hvug). However, not only such solitons are essential at low tem- 
peratures but also they are known as the only static solitons allowed at equilibrium. 
At low temperatures, solitons' widths increase as a function of bias voltage V and 
temperature T (an outcome of a decrease of solitons' amplitude at low temperatures 
as discussed in ref. [5D). Detailed behavior of the self-consistent equation (|4-4p and 
current between two reservoirs (j4-5p were particularly studied in ref. |8|). Since the 
amplitude of solitons follows the same self-consistent equation as the isotropic lattice 
distortion, the stabilities studied in ref. [5]) can be considered to be the same as the 
stabilities of the amplitude degree of freedom of the solitons. 

Solitons' amplitudes or widths {A{x) = hvug) characterize fermionic transport 
because current between reservoirs only depends on their amplitude and width. To 
be more precise, higher bias voltage V and/or higher temperature T imply smaller 
amplitude (larger width soliton) and lower current at low temperatures. Detailed 
relationships between temperature, bias voltage and solitons' amplitude can be un- 
derstood by replacing isotropic lattice distortion with solitons' amplitude, as in the 
previous paper. The results of the isotropic lattice distortion, i.e., amplitude of the 
solitons, are roughly summarized as follows. 
• at constant bias voltage 

The phase transition between ordered and normal phases could be first or 
second order depending on the parameters. In some parameter regions, the 
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voltage-current characteristics are S-shaped (namely, there is negative differen- 
tial conductivity). Negative differential conductivity appears when temperature 
is lower than a certain threshold value. 
• at constant current 

All the non-trivial solutions of the self-consistent equation are stable and the 
phase transition between the ordered and normal phases is always second order. 
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Appendix A 

Derivation of a continuous model 

In this appendix, we explicitly show the Green function in the soliton case. For 
the TLM model, the Green function obeys 



D+g+a{x, y, u) = ujg~a{x, y, uj) - 6^aSix - y) 
D^g^^{x, y, uj) = ug+^{x, y, uj) - 6+aS{x - y) , 

where D± is defined by D± = + A{x)/h. Then, g-^{x,y : u) follows 



D_^D^g_^{x,y : uj) 



^2 ^ V dAjx) ^ A{x)^ 
dx"^ h dx h? 



g-a{x,y : w) , 



at X ^ y. For the soliton A{x) = AQtaiihKa{x — a), the operator Dj^D^ is easily 
caluculated; 



q2 ^ ^2 

D^D^ = —v'^——;t H — —VKs sech^K^fx — a) + — tanh^sfa; — a) 
ox^ n h'^ 



In general, the solutions of (jA-ip are expressed in terms of the hypergeometric func- 
tion. In this paper, we restrict ourselves to the case of Aq = hvus- As such, the 
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green function is given by 

' {{hvKf + A{0)A{x)} sinK(x) + hvK{A{x) - A{0)} cosk{x) 



{hvYtiujDiuj) 

A{y) sin — I) — hvK cos K{y — I) {x < y) 

A[x) sin k{x — i) — hvK cos k{x — i) 
{hv)'^KUjD{u>) 

X {{hKv)'^ + A{0)A{y)} sin Ky + hvK{A{y) - A{0)} {x > y) 
u> [hvK COS nx + A{Q) sin nx) sin K{y — 



g—ix,y;u;) 



v'^kD{uj) 

CO sin k{x — £) [hvK cos ny + A{Q) sin Ky) 



v'^nD{uj) 



{x < y) 
{x > y) 



g+^{x,y;u) = g^+{y,x;uj) 



\^{hvK)^ + A{x)A{0)^ sinKx + hvK\^A{x) — A{0)^ coskx sin K{y 



hv'^KD{(jj) 

(Z\(x) sin k{x — i) — hvK cos k{x — ^)) (^(0) sin Ky + hvK cos ny^ 

hv'^KD{io) 

where k = \J (huoY — A'q/ (hv) and D{uj) = hvncos + A{Q) sinK^. 



{x < y) 
{x > y) 
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